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1 Introduction 



The search for a mechanism explaining electroweak symmetry breaking has been a major goal for many 
years, in particular the search for a standard model (SM) Higgs boson, see for instance Ref. [T] and Ref. [2J. 
As a result of this an intense effort in the theoretical community has been made to produce the most accurate 
NLO and NNLO predictions, see Refs. [3141516] . There is, however, a point that has been ignored in all these 
calculations: the Higgs boson is an unstable particle and should be removed from the in/out bases in the 
Hilbert space, without destroying the unitarity of the theory. Therefore, concepts as the production of an 
unstable particle or its partial decay widths do not have a precise meaning and should be replaced by a 
conventionalized definition which respects first principles of quantum field theory (QFT). 

The quest for a proper treatment of a QFT of unstable particles dates back to the sixties and to the work 
of Veltman [8] (for earlier attempts see Ref. [9]); more recently the question has been readdressed by Sirlin 
and collaborators [10] . Alternative approaches, within the framework of an effective theory can be found in 
Ref. [H]. 

In this paper we discuss the relation between physical observables and Higgs pseudo-observables by 
considering the extension of their definition to higher orders in perturbation theory, confronting the problems 
that arise when perturbative corrections in quantum electrodynamics (QED) and quantum chromodynamics 
(QCD) are included. Numerical results are also presented. Our work can be seen as an extension of 
complex-mass schemes to include complex external momenta (for previous work see also Ref. [12]). addressing 
systematically the question of the analytical continuation of Feynman loop integrals. 

This paper is organized as follows. In Section [5J we summarize the conceptual setup. In Section [3] we 
present general arguments on complex poles. In Section [4] and [5] we discuss pseudo-observables, on-shell 
observables and unitarity. The analytical continuation of Feynman loop integrals into the second Riemann 
sheet of the S -matrix is examined in Section |6] In Section [7] we present the inclusion of QED and QCD 
corrections and renormalization schemes are highlighted in Section[8l Numerical results are given in Section[9] 
and in Section [TU] we close with our conclusions. 

2 Formulation of the problem 

There are two old questions in relating measurements to theoretical predictions: 

- Experimenters (should) extract so-called realistic observables from raw data, e.g. o~(pp — ¥ 77 + X) and 
need to present results in a form that can be useful for comparing them with theoretical predictions, 
i.e. the results should be transformed into pseudo-observables; during the deconvolution procedure one 
should also account for the interference background - signal; 

- Theorists (should) compute pseudo-observables using the best available technology and satisfying a list 
of demands from the self-consistency of the underlying theory [14] . 

Almost from the start it is clear that a common language must be established in order to avoid misun- 
derstandings and confusion. A typical example can be found in Higgs physics where, frequently, one talks 
about Higgs production cross section or Higgs partial decay widths. After the discovery phase, in absence 
of which the future of high energy physics cannot be ascertained, one will need to probe the properties of 
the discovered resonance, like spin and couplings. In this case different sources will start talking about the 
same thing but with different languages. We will indicate a reasonable language within the context of a 
perturbative expansion of a gauge-invariant QFT in this paper. 

The Higgs boson, as well as the W or Z bosons, are unstable particles; as such they should be removed 
from in/out bases in the Hilbert space, without changing the unitarity of the theory. As mentioned before, 
concepts as the production of an unstable particle or its partial decay widths, not having a precise meaning, 
are only an approximation of a more complete description. The inconsistencies associated with the on-shell 
LSZ formulation of an unstable external particles become particularly severe starting from two-loops, as 
described in Ref. [5]. 
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Suppose that we want to combine a Higgs production mechanism, say gluon-gluon fusion, with the 
subsequent decay H — > 77. The process to be considered is, therefore, pp — > 77 + X and it is made of a part 
that defines the signal, e.g. 

PP -> 99(^ H -> 77) + X, (1) 

and by a non-resonant background. The question is: how to extract from the data, without ambiguities, a 
pseudo-observable to be termed Higgs partial decay width into two photons which, at the same time, does 
not violate first principles? Once again, one should be aware that there is no Higgs boson in the in-state, 
therefore the matrix element < 77 out|i? in > is not definable in QFT and this ill-defined quantity should 
be replaced by a pseudo-observable which closely resembles the intuitive concept of a decay width, can be 
unambiguously extracted from the data and respects all fundamental properties of the theory; in this way 
we replace a non existing observable with a conventional definition. A proposal in this direction can be 
found in Ref.[10j; here we revise the proposal, improving it by considering the extension to higher orders in 
perturbation theory, confronting the problems that arise when QED(QCD) corrections have to be included 
and present numerical results for Higgs physics. 

At the parton level the S -matrix for the process i — > f can be written as 

Sfi = Vi(s) A H (s) Vf(s) + Bif(s), (2) 

where Vi is the production vertex i — > H (e.g. gg — > H), Vf is the decay vertex H — > f (e.g. H — > 77), 
A H is the Dyson re-summed Higgs propagator and Bif is the non-resonant background (e.g. gg — > 77 
boxes). In the next section we will introduce the notion of complex pole. A vertex is defined by the following 
decomposition [16] , 

Vf(s) = Y,V f a (s, {S}) Ff({p f }) (3) 

a 

where s = —P% (with P H = YlfPf)) s © {&} 1S the set of Mandelstam invariants that characterize the 
process H — > f, V? are scalar form factors and the Ff contain spinors, polarization vectors, etc. 



3 The complex pole 



In this section we introduce the notion of the complex pole [15] following closely the analysis of Ref. [17] . 
Let Ai be the lowest order propagator for particle i and A.; the corresponding dressed propagator, i.e. 



Ai = - 



(4) 



Let us analyze in more details the definition of the dressed propagator: to begin with, consider a skeleton 
expansion of the self-energy S = 16 7r 4 i S with propagators that are resummed up to O (n) and define 



A^(s) = Af 00 [l + Af( S ) 4 l) (• , Af" ^OO) 



(5) 



where, omitting an overall factor — z/(2 7r) , the Born propagator (tensor structures are easily included) is 



A, (0) ( S ) 



1 



,2 ' 



(G) 



If it exists, we define a dressed propagator as the formal limit [5] 



A i (a)= lim AHs), Aj(s) = - A- j (s) 1 + A u ^( s ) E« (s , A,(s)) 

which coincides with the Schwinger-Dyson solution for the propagator. 
The Higgs boson complex pole (s H ) is the solution of the equation 



(7) 



s H - Ml + Shh(s h ) = 0, 



(8) 
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where M\ is the renormalized Higgs boson mass; here we assume that all counter-terms have been introduced 
to make the off-shell Green's function ultraviolet finite, respecting locality of the counter-terms. We now 
examine more carefully the self-energy to all orders in perturbation theory since, often, there is some confusion 
with statements that are formulated to all orders and applied to a truncated perturbative expansion. Now 
consider the, all orders, self-energy, 



X HH (s,M 2 H ,0 = 54"2(s,M2,£) ff 2 ", 



(9) 



where £ is the gauge parameter (extension to more than one gauge parameters is straightforward) and g is 
the renormalized coupling constant. From arguments based on Nielsen identities, see Ref [16] . we know that 
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s H = 0, 



(10) 



i.e. the location of the complex pole is £ independent; as a consequence the self-energy is £ independent 
too, since the two differ by a renormalized quantity, obviously £ independent. We consider first the one-loop 
approximation for the self-energy: from its explicit expression we are able to derive the following relation: 



(11) 



E« (a, M 2 B , = Ei 1 ! . tl {s,M*) + (8- M 2 H ) <S> H (s, Mf , £)■ 
where, in a general gauge, we obtain 

Ml $ H = -1 1 {s+Ml) [B d (s, M 2 z , M 2 z ; 6) + 2 B d (s, M W ,M W ; U)] +2A d (M z ^ z )+4 A d (M w , f w ) 

B d (s,m,m,£) = B (s,£m,£m) - B (s, m,m) , A d (m, £) = A (£ m) - A (m). (12) 

The symbols Aq,Bq, etc. are the usual scalar, one- loop functions. 

It needs to be stressed that the splitting between gauge dependent and gauge independent quantities 
is only defined modulo a £ -independent constant. Our definition of the invariant part is that it coincides 
with the expression in the 't Hooft-Feynman gauge (i.e. £ = 1). Furthermore, finite renormalization (i.e. 
replacing renormalized parameters with a set of experimental data points after having removed ultraviolet 
poles by means of local counter-terms) amounts to replace 



showing that 



Ml = s H + X HH {s H ,Ml,0, 



7T7 ^hh(Sh, ShtO — 0' 

<9£ 



(13) 



(14) 



so that, at one-loop, the Higgs complex pole is gauge parameter independent if the self-energy is computed 
at Ml = s H , the basis of the so-called complex-mass scheme (see Ref. [TB] and also Ref. [T7]). From Eq. ([TU)) 
and from the one- loop result in Eg. dill) , we derive the following 



EL"2 (s h , Ml, = E w . , (s H , Ml ) + E w . 4 (s„ , Ml , £) , 
X%l ;S (s s ,Ml,Z) = Z^}(s H ,Ml)$ s (s s ,Ml,t;). 

Using Eq. (IT5l) we can rewrite Eq.® at the two-loop level as 

1V1 H — s h ~r 9 ?->hh ;i{S H , S H ) + h ; i\ s h t S H ) ~t Z^ HH . j(S H , S H ) 



dMl 



E H /f ; j (S H , M H ) 



M%=S H 



(15) 



(16) 



We find as worthy note that Eq. (|16[) can be easily generalized to all orders in perturbation theory showing 
that, order-by-order, the gauge dependent part of the self-energy drops out in the equation for the complex 
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pole of the particle. The complex pole, sitting on the second Riemann sheet of the S -matrix, is usually 
parametrized as 

It is worth noting that a consistent treatment of external (s) and internal (M%) masses allows the extension 
of the complex mass scheme beyond one- loop, without the need of expanding the self-energy around s H = fi%, 
as frequently done in the literature. In partial contrast to the traditional complex mass scheme, Ref. j!8j . in 
our approach (described in Ref. [17] ) it is the finite renormalization equation and not the Lagrangian that 
is modified. Indeed, calling the scheme complex mass scheme is somehow misleading; to the requested order 
we replace everywhere the renormalized mass M\ with s B + Y} BB (s B ) which is real by construction; if only 
one-loop is needed then Mj — > s B everywhere, therefore justifying the name complex mass. 

The quest for gauge invariance and the consequent introduction of a complex pole instead of an on-shell 
mass signal has a certain degree of ambiguity in defining the Higgs boson mass (as well as the mass of any 
unstable particle). The most convenient choice, for all practical purposes, is represented by the square root 
of the real part of s H , although 

& = 04 +7!) 1/2 (is) 

also has several advantages [H] and will be used in our numerical results. 

There is a final comment for this section: the complex pole for an unstable particle, parametrized 
according to Eq.®, must correspond to a negative imaginary part; otherwise, even the Wick rotation 
cannot be safely performed. Consider the case of the Higgs boson, ii channels that do not satisfy the 
negativity condition for the imaginary part below the 4 mf (real) threshold are excluded in the evaluation 
of s H . As we already mentioned the contribution to the imaginary part of s H from a given channel below 
the corresponding real threshold (WW, ZZ and tt) represents an approximation to the corresponding 4/ 
and 6/ cuts, i.e. H — > WW, ZZ —> 4/ etc, which is acceptable only when the corresponding 7^ is positive, 
a condition which fails at one-loop for tt intermediate states when the top quark mass is kept real; in this 
case tt intermediate states never contribute, in our scheme, to j H below threshold, i.e. they are discarded. 
It is interesting to note that this problem completely disappears if we allow for a top quark complex pole 
(instead of real on-shell mass). Numerical examples will be discussed in Section |9j unfortunately the top 
quark total (on-shell) width is poorly known, therefore inducing large uncertainties on the corrections. In 
the numerical analysis we use T t < 13.1 GeV, based on the experimental upper limit of Ref. [20] . 



4 Extracting a partial decay width 

In this section we examine our options to define a pseudo-observable which is related, as closely as 
possible, to a realistic cross section and shares as many features as possible with the corresponding on-shell 
definition of a partial decay width. If we insist that \H > is an asymptotic state in the Hilbert space then 
the observable to consider will be < / out | H in > , otherwise one should realize that for stable particles the 
proof of the LSZ reduction formulas depends on the existence of asymptotic states 

I pin >= lim [ d 3 xH(x)iB t e lpx \0 >, (19) 
t— >— 00 J 

(in the weak operator sense). For unstable particles the energy is complex so that this limit either diverges 
or vanishes. Although a modification of the LSZ reduction formulas has been proposed long ago for unstable 
particles, see Ref. [21) . we prefer an alternative approach where one considers extracting information on the 
Higgs boson directly from 

< / out I H >< H I i in > + V" < / out | n >< n \ i in >, (20) 

for some initial state i and some final state / and where {n} © H is a complete set of states (not as in the 
in/out bases). As we are about to see, the price to be paid is the necessity of moving into the complex plane. 
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Define H HH (s) as 



IWs) = 



then the, Dyson re-summed, Higgs propagator becomes 

A HH (s) = (s - sx)- 1 l + n HH (s) 
Using Eg. (|2"21 we can write Eq.J2]) as 



z H = i + iw 



Sfi - 



Z- 1/2 (s)V f (s) +B if (s). 



From the S -matrix element for a physical process i — > f we extract the relevant pseudo-observable, 

S(H C ^ f) = Z- 1/2 ( SH )V f (s H ), 
which is gauge parameter independent - by construction - and satisfies the relation 



•>fi 



S (i -> H c ) S (H c -> /) 



non resonant terms. 



The partial decay width is further defined as 

(2tt^ 4 



fi a r (H c -+f)= J d* f (P H , { Pf }) J2 \s (H c -> /) 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



spins 



where the integration is over the phase space spanned by \f >, with the constraint P H =J2 pf. One should 
not confuse phase space and the real value of s = — P|, where the realistic observable is measured, with the 
complex value for s, where gauge invariant loop corrections must be computed. The choice of (phase 
space) where to define the pseudo-observable is conventional, e.g. one can use the real part of s H . Indeed, 
the r.h.s. of Eq. (l24|) satisfies the property 



d_ 



Z- 1/2 ( SH )V f (s H )=0 



to all orders in perturbation theory. If we define 

oo 

V f (s,Ml)=J29 2n+1 



vj:](s,Ml)+V}%( S ,Ml) 



(n) 



(27) 



(28) 



we obtain, expanding in powers of the coupling constant <?, that 
vgl(a H ,a B ) = ±vW$ H (s H ,s H ), ^(s„,s„) = -l$ H (s H ,s H )\vjll(s H ,s H )-~vW$ H (s H ,s H )\, (29) 



etc. One last example of a basic fact: Nielsen identities give the structure of the gauge parameter dependent 
vertex and self-energy order-by-order in perturbation theory. It is important to stress at this point that the 
renormalized mass should be replaced consistently with the use of Eq. (|13|) . 

To summarize, only s„ is a meaningful quantity and a definition of the real mass or of the total width is 
conventional. From Eq.© one has 

HhJh = ImE ffiJ (%), (30) 
and it should be evident, from Eq. (|26|) . that 7„ ^ J2f r {H c /)• The reason can be understood when 
we consider a simple example, a toy model with C- ln t = m? </)<t + er~ (with massless a -particles). Already at 
one-loop, we find 



ImE^(s) 



16 7T 



■7T, ImE^s^) 



16 



i 2 f 1 



arctan ( — 
M0 



(31) 



5 



While the first relation in Eq. (j3"T1) (real s) satisfies the cutting equation [7] the second (complex s) does not. 
For a proper perspective it must be recalled that when we expand, 'E HH (s H ) = E flH (^) + . . . , the cutting 
equation is restored at NLO but it will still be violated at NNLO, as pointed out in Ref. [IB]- Therefore, our 
conventional definition of the Higgs total decay width will be r tot (i? c ) = J^f T (H c —> f). 

To set the stage, it may be well to recall that the breakdown of a process into products of pseudo- 
observables can be generalized to include unstable particles in the final state; an example is given in Fig. Q] 
where the (triply-resonant) signal in gg — > 4f is split into a chain gg —> H (production), H — > W + W~ 
(decay) and W — > ff (decays). 




S w 

Figure 1: Gauge-invariant breakdown of the triply-resonant gg — > 4f signal into gg — » H production, H — > W + W 
decay and subsequent W —¥ ff decays. 



5 Pseudo-observables, on-shell observables and unitarity 

When we consider all the possible decay channels of an on-shell standard model Higgs boson we obtain 
that up to an on-shell mass m H « 140 GeV the Higgs boson is very narrow while the width rapidly increases 
after the opening of the WW and ZZ channels. 

Even this statement should be carefully examined since W and Z bosons are unstable particles to be 
removed from the in/out bases of the Hilbert space. For real W, Z masses the Higgs boson width is related to 
the cuts of the self-energy and the statement under examination is based on the (say one-loop) two-fermion 
cut, two-boson cut, etc. 

Unitarity follows if we add all possible ways in which a diagram with given topology can be cut in two 
separating S from S>. For a stable particle the cut line, proportional to the positive energy part of the 
propagator, contains a pole term 2 i ir 6(po) S(p 2 + m 2 ), whereas there is no such contribution for an unstable 
particle. We express Im £ in terms of cut self-energy diagrams and repeat the procedure ad libidum, therefore 
proving that cut unstable lines are left with no contribution, i.e. unstable particles contribute to the unitarity 
of the S— matrix via their stable decay products [5]. 

From this point of view the second cut of the Higgs self-energy (after the two-fermion cut) is the four- 
fermion cut, not the two-boson one (once again, the cutting of a line corresponding to an unstable particle 
contains no pole term). How bad is the choice of cutting two, stable, W boson lines with respect to cutting 
four fcrmion lines and summing over all fermions, i.e. how bad is the on-shell approach, at least from a 
numerical point of view? 

If one evaluates the ratio 

r (H -> VV) BR (V ->■ 2f) BR (V -> 2/') / T (H 2/ + 2/') (32) 

the results of Ref. 22 show that the on-shell phase space for the WW or ZZ final state introduces an error 
of the order of 10% near the threshold, which is still satisfactory. Using the complex mass scheme which, 
in turns violates unitarity, will improve upon the on-shell result since the internal V masses are themselves 
complex poles. Remarkably, the complex mass scheme represents a method which is, at the same time, 
predictive and gives the best available approximation to the use of a full (Schwinger-Dyson) re-summed 
theory, a formal solution of the problem which, however, poses an insurmountable barrier for the technology 
of today. 
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6 Loop integrals with complex masses and invariants 

In this section we analyze the correct definition of Feynman integrals with complex masses and Mandel- 
stam invariants. On a more formal bases one should say that unstable states lie in a natural extension of the 
usual Hilbert space that corresponds to the second sheet of the S -matrix; these states have zero norm and, 
therefore, escape the usual prohibition of having an hermitian Hamiltonian with complex energy |21) . On a 
more pragmatic level we use the guiding principle that Green's functions involving unstable particles should 
smoothly approach the value for stable ones (the usual Feynman — i prescription) when the couplings of 
the theory tend to zero. 

The whole problem can be summarized as follows: in the limit of zero couplings all particles are stable 
and we define Green's functions in the cut s -plane, where s is the selected invariant to be continued into the 
complex plane. For the free theory of stable particles, according to Feynman prescription, the value of the 
argument of some function lies, say below the cut (which coincides for example to the real negative axis); 
during continuation of s we may cross the cut, which means that we have to continue the function into its 
second branch. 

For the simple case that we have just described the Green's function is then defined through its value on 
the principal branch in all quadrants but the second, where continuation to the second branch is required. 
This new function will have a cut on the positive imaginary axis and special problems may occur, especially 
when we want to do analytical continuation at the level of integrands and also internal masses in a given 
Feynman diagram are complex, as required by any realistic complex-mass scheme. Green's functions are 
given in terms of Feynman parametric integrals and our main point will be: how to define the same integrals 
but properly continued to complex internal masses and complex external invariants? One of the difficulties 
of the problem lies in having masses and invariants complex at the same time which introduces subtleties in 
the analytical continuation which are not present if, say only masses or only invariants are made complex. 

6.1 General setup 

To start our analysis, consider a scalar <fia 2 theory with > 2m CT , i.e. <f> is unstable; the <j> propagator 
(with s = — p 2 ) is 

(33) 

where factors (2 tt) i have been omitted. The inverse function, A _1 (s) is analytic in the entire s -plane except 
for a cut from s — 4 m 2 to infinity along the real axis. The function is defined above the cut, A -1 (s +i0) 
and the analytical continuation downwards is to the second Riemann sheet, i.e. 

Aj^s-iO) = A -1 (s + iO) = A- 1 (s-iO) + 2iTrp(s), (34) 

where 2 in p(s) is the discontinuity of the function across the cut. For a complete discussion of the analytical 
continuation see, e.g., Ref. [2"4"] . 

We need a few definitions which will help the understanding of the procedure for the analytical continu- 
ation of functions defined through a parametric integral representation. The logarithm is defined by 

ln (fe) z = ln<- 0) z + 2ink, k = , ±1 , ±2 , . . . (35) 

where In' ' z denotes the principal branch (—tt < arg(z) < +7r). From now on we will omit the superscript 
that denotes the principal branch of the logarithm. Let z± = zq ± i and z — z R + i Zj, we define 

in [z,z ± )-\ \nz±2iTr9(~z n ) e(TZj), W 

i.e. the first Riemann sheet for all quadrants but the second where the function is defined in the second 
Riemann sheet. 

Our first definition of the In^ -functions in Eq.(|36j) is the most natural in defining analytical continuation 
of Feynman integrals with a smooth limit into the theory of stable particles; the reason is simple, in case 



A = 



a - M 2 + £^( S ) 
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some of the particles are taken to be unstable we have to perform analytical continuation only when the 
corresponding Feynman diagram, in the limit of all (internal) stable particles, develops an imaginary part 
(e.g. above some normal threshold). However, in all cases where the analytical expression for the diagram 
is known, one can easily see that the result does not change when replacing zq with z R , the second variant 
in Eq.(|36]). 

As we are going to discuss in the following sections there are cases where one would like to perform 
an analytical continuation at the level of integrand in the Feynman parametric representation of a given 
diagram; often the integration contour has to be distorted into the complex plane with the consequence that 
z R 7^ zq and sign(z Ji ) ^ sign(zo)- In this case we need a more general definition of Ivi^: 

Definition: Let z(T) £ C (T £ R) be an arbitrary complex function of T; when we want to continue 
zo = z(0) (not in the second quadrant) to Zf — z(Tf) we must look for a real T c with < r c < Tf 
such that z c — z(T c ) is real and negative (for simplicity we assume the case a monotonic zr = z(T)): 
then, VT : r > r c we replace lnz with its analytical continuation into the second Riemann sheet, 

ln ± (z r ; z ) = In z T ± 2 i ir 6 (-Re z T ) 9 (=Flm z r ) • (37) 

For all practical purposes Eq. (|37|) can be replaced with the second variant of Eq. (f36| (with z — > zp) 
which, from now on, will be our definition of analytical continuation. 



6.2 Analytical continuation of the Euler dilogarithm 

We consider now the Euler's dilogarithm, Li2 (z); if we denote its principal branch by hi^' \z) (0 < 
arg(z — 1) < 271"), than for any branch (see, e.g. [25]) we have 

U^" m) (z) = U { °' 0) (z) + 2 nn i ln {0) z + imn 2 , n , m = , ±1 , ±2 , . . . (38) 
The question that we want to analyze is the following: given 

f 1 dx 

Li 2 (M 2 +i0) = - / — ln(l- M 2 x-i0) , ImLi 2 (M 2 + i 0) = tt In M 2 6 (M 2 - l) , (39) 

how do we understand Eq. (|38|) in terms of an integral representation? Let us consider the analytical contin- 
uation from z + = M 2 + i to z = M 2 — iMT and define 

1 dx 

— In" (l - zx; 1 - z + x) . (40) 
x 

With x(x) = l-zx = l- (M 2 - i MT) x, we have x(0) = 1 and = (l - M 2 , M V) . If M 2 > 1 we have 
that x crosses the positive imaginary axis with Im% = T/M. As a result we obtain 

I = L4 ' 0) (z) + 2in In M 2 , (41) 



which is not the expected result since I does not reproduce the correct continuation of Li 2 given in Eq. (|38|) . 
The mismatch can be understood by observing that ln~ has a cut along the positive imaginary axis (of x) 
and, in the process of continuation, with x £ [0, 1], we have been crossing the cut. Nevertheless, we insist on 
defining analytical continuation at the level of integrand, instead of working directly on the result, because 
it is the only practical way of dealing with multi-loop diagrams where an exact result is not known. The 
solution consists in deforming the integration contour, therefore defining a new integral, 

I G = f —kT(l-zx;l-z + x), (42) 
J c x 

where the curve C is given by two straight segments, < x < 1/M 2 — e and 1/M 2 + e<£<l(e— ^ + ), 
plus a curve C defined by 

C\u) : {x^u + i 1 -^}, J^<u<^ (43) 
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The integral over C" is downwards on the first quadrant an upwards on the second (along the cut of In ) . 
Integration of In - over C gives —2 i ir (In M 2 — In z), showing that 

L#< 0) (z)=I o , (44) 

the correct analytical continuation. Therefore we can extend our integral, by modifying the contour of 
integration, to reproduce the right analytical continuation of the dilogarithm. 

6.3 Continuation of analytical results 

Having introduced a simple example, we consider now one-loop two-point functions where both masses 
and the external invariant are made complex. Let 

x{ x ) = s p x2 + ( m 2 — m i ~ s p) x + to ij (45) 

s P = M 2 —iTM, m 2 = pi 2 - i ji m. (46) 
The function B$ is originally defined, for real s P and real (equal for simplicity) internal masses, by 



i r 1 

B [M ; fj,,{j) = - - / dx ln(x - i 0), (47) 
o 



where e" 1 = 2/(4 — n) — ^ E — \mr (j E w 0.5772 being the Euler-Mascheroni constant), and we need the 
analytical continuation to arbitrary values of s P (i.e. M 2 — > M 2 — iMT with T > 0); we assume, for a 
moment, real internal masses (7 = 0) and M 2 > 4/z 2 ; the analytical result is 

B (M 2 ;n,v) = T-ln4 + 2-/3 1n§±± (48) 

where fi R is the renormalization scale and (3 2 = 1 — 4/x 2 /M 2 . For the continuation to M 2 — > M 2 — iMT, 
fj 2 —> /i 2 - i 7 /i we have to compute the logarithm of z UST = z^ ST + iz^ ST , which is a function of T,j 
(interacting theory of unstable particles) . Let 

4 T = lim zJf T ±z0 = 4 T ±*0, (49) 
r,7— >o 

where the ±i follows from Feynman prescription /1 2 —> /i 2 — i 0. We use the second variant of Eq. (|36| and 
define 

ln ± (z UST ; z^ T ) = In z UST ± 2 i n 9 {-z v R ST ) 9 (t< st ) , (50) 

which satisfies 

lim In* (z UST ; zf) =lnzf, (51) 

and it is equivalent to have lnz on the second Riemann sheet, but only when z is continued into the second 
quadrant. There is one awkward possibility; it corresponds to starting from zi T with z^ T < and requiring 
continuation to z" ST > and z^ ST > 0. Using Eq. (|48|) we derive 

r2. ... _^ ,1 ^ , o a ,_- [Pc + l . 13 + 1 



B (M 2 ;m,m)^ T -ln— +2-/3 c ln- ^— ; ^— (52) 

where /3 2 = 1 — 4 m 2 / s, p . It is worth noting that there is never a problem when internal masses are real and we 
continue to complex p 2 . Otherwise we first continue to complex internal masses using the fact that internal 
(complex) squared masses have a negative imaginary part. Consider this continuation for one-loop diagrams: 
with L -external legs we can always fix a parametrization where the coefficient of m 2 is 1 — x\, the one of 
mf is Xi—i — Xi, up to m 2 which has coefficient x L -\ where the parameters satisfy < i L _i < ... < x\ < 1 
(i.e. all coefficients are non- negative) . Less straightforwardly the same holds for multi-loop diagrams. 
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Then we continue from T = 0; considering Eq. (j52")) and denoting by £ the ratio (f3 c + l)/(j3 c — 1) we have 
Re/3;? > (in the region above real thresholds) and Im/3 2 > for T = 0. At T — (M/fi)j /3 2 crosses the 
positive real axis from above; this corresponds to Q crossing the cut where we move into the second Riemann 
sheet of the logarithm of Eq. (j5"2")) . After that one has Im£ > and the forbidden region is reached when 
Re£ > 0, which corresponds to | f3 c \> 1. Once again, for 7 = the latter is never satisfied. In general, the 
condition requires solving a cubic equation in T with only one real, negative, solution. The forbidden region 
requires, therefore, T < 0. 

To continue our analysis of one-loop functions, where analytical results are known, we only need to define 



Li J u 



' 5T ; z| T ) = Li 2 0) T 2 i ir 9 (zjf T -1)6* (±z, UST ) In z l 



(53) 



For our purposes, namely for the processes that we are considering, we only need one additional function. 
The most general scalar three-point function that is needed will be 



Co (0,0,p 2 ; mi ,m 2 ,m 3 ) = ^{E 



with four different roots 



x = l 



i=l,3 



£3 



Li 2 



Xq 



1 



X - Xi 



-Li 2 



Xq - Xl 



+ lnx r\ (x\ —xo,%2- £0)}, 

(54) 



in-. 



p 2 + m\ — m 2 =F A 1//2 (—p 2 , m i, m?) 
xi,2 = ^ , (55) 



where A is the Kallen function. Analytical continuation requires the replacement Li 2 — > Li 2 with limiting 
(free theory of stable particles) cases given by 



x\ — > x\ — i 0, 



X2 -» X2 + iO, 



£3 — > 2^3 — i sign(m 2 — m\) 0. 



(56) 



As a final observation, there is no need to continue the square root /3 C in Eg. (j5"2"j) below threshold (f3 2 < 0) 
since in this case j3 c is imaginary and the change of sign when we move from the principal root is compensated 
in the product /3 C times the logarithm. 

Finally, in Eq. (l54"|) and for one- loop processes with more scales and more than three legs one has to 
introduce a generalization of 't Hooft-Veltman 7/ -functions [29 on the second Riemann sheet. The definition 
is as follows: 

In - (x y) = In - x + In - y + rf~ (x, y) , (57) 
ff{x,y) =2nr{6(x 1 ) [e(-x R ) - %,) 0(-^)] + 6{ Vl ) 6{-y R ) + 9{z t ) [0(-x T ) 9(- yi ) - 0(-z R )] }, (58) 
with z = xy and z — z R + i z, etc. 



6.4 Continuation at the integrand level 

We now turn to analytical continuation at the integrand level, according to our procedure where all 
Feynman integrals are treated according to their parametric integral representation. 

Let us consider the specific example of the previous section: a scalar two-point function corresponding to 
two internal equal masses, m 2 = /i 2 — i /j,-/ and s p — M 2 — i AIT. Due to the symmetry of x( x )> m Eq. (|47j) 
the integral with < x < 1 can be written as twice the same integral with < x < 1/2; the argument of the 
logarithm goes from Rex = M 2 > to Re% = /i 2 — M 2 /4 < (above threshold) with Imx = — iO. We have 
to define the analytical continuation M 2 — > M 2 — i MY; since, for any x, x cannot cross the cut, it must be 
analytically continued into a second Riemann sheet above the cut. A similar situation occurs for complex 
internal masses: the integration with respect to x is 

x = (/A -ifJL-y) x= (V- \m\ -i^+'-MvY (59) 
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Let X = x (1 — x) with < X < 1/4, select a value for x, when T > (fij)(MX) continuation is into the 
second Riemann sheet. Of course, for M 2 < 4/j, 2 , \ remains on the first Riemann sheet for all values of F. 
The variable x is such that 



Rex = 
Imx = 



for x — R± = — 



for x = 1+ = — 



l±\ 1 



M 2 



l±Wl-4 



MT 



(60) 



The second equation requires MT > 4 ^7 for I± to be real and £ [0,1]. At x — I± the condition Re x < 
requires T/x < M7. Therefore, for those values of T and x that satisfy the conditions 



4-^7<T< — 7, 



J_ < x < I 4 



R-<x< R 4 



we have Rex < 0, Imx ^ an d hix must be continued into the second Riemann sheet. 
The new definition of the Bo -function is as follows: 

1. 1 



^ In (x;x_), 



X- = X 



r,7i=o 



iO, 



(61) 



(62) 



Different possibilities are illustrated in Fig. [2] where we show x( x ) f° r two equal (complex) internal masses. 

Im \ ' 




Figure 2: Analytical continuation from real p 2 to complex p 2 as seen in the x-pl ane with x{ x ) ~ — s p x(l — x) + 
H 2 — ifJ,J, with s P = M 2 — iMY and x £ [0, 1]. Solid lines represent the continuation for a low value of M with a 
very small value for T. With increasing values for M we reach the situation illustrated by the dot-lines, \ moving 
into the second quadrant, i.e. \ on the second Riemann sheet. Case 1 holds for V < (M//i) 7 whereas case 2 holds 
for r > (M/fj,) 7. Black circles correspond to x = 0, x — 1 whereas white circles correspond to x — 1/2. 

In any realistic application the complex pole equation returns, for low values of M, small values of T and Re x 
is always positive, never requiring analytical continuation into another sheet; when M increases T increases 
too and we find values of x that requires the continuation In — > ln~ . This will happen for x > 1- in case 1 
(which requires T < (M/fijj) and for x > R_ in case 2 (which requires T > (M//j,)j). 

The same example can be discussed in the x complex plane; in this case, when M > 2 /i and T = 7 = 0, 
the cut is on the real axis between R_ and R + (Eq. ([6H)) ) and the integration is < x < 1/2. The integral is 
originally defined (T, 7 = 0) above the cut, i.e. for x + i 0. Analytical continuation means that for increasing 
imaginary parts we reach a point where the integration path is continued into the second Riemann sheet (at 
T = [Mj fi) 7, we have /_ = i?_ and, for higher values of T, the continuation to the second Riemann sheet is 
required as soon as the cut is reached). From this point of view the integral is better understood in terms of 
a variable z{x) = u + iv, such that x = —M 2 z (1 — z) + fi 2 and the integration is performed along the curve 



u(l- 



-2U) = —jr. 

; M 2 ' 



Z = fj,j-MTx(l-x), 



1-U 



U 4 



Ax{l-x)-l\U 2 -A Wi =Q. (63) 



11 



Note that x = I_ corresponds to z = real. For real internal masses we have z(0) = and z(l/2) = 
(r/(2M)) 1 / 2 (1 —i/2). In the z -plane the logarithm has a cut on the positive real axis between R_ and 

It is worth mentioning that case 2 of Fig. [5] corresponds to an integration path that crosses the cut of 
In - across the positive imaginary a; -axis, similar to the case of the dilogarithm discussed above. Therefore, 
the correct analytical continuation, for case 2, goes as follows: the integration path in z -space (Eq. (|63)0 is 
moved into the complex plane and goes into the lower half-plane instead of reaching the cut of the logarithm 
(which is between i?_ and R + , see also Fig. |3]). 

In order to insure that the analytically continued integral has a smooth limit r, 7 -> we deform the 
integration path by insisting that the cut (of In) must be crossed at z = R- (note that for case 2 we have 
J_ < i?,_) where we perform a continuation into the second Riemann sheet. In this way we add to Bq (on 
top of a factor — 2iirf3) a new contribution which is easily computed in the x -plane and it is related to the 
discontinuity of In \ along a curve C parametrized by 



C(t) : {*=!-*+»/(*)}, 



2 I M 



/(*)=o -TT?* + 



1 



1/2 



}■ 



(64) 



with (3 2 = 1 — 4/i 2 /M 2 > and where (3 < t < (3; here (3 is the value of t where Rex(i) = ImxW = 0. 
The integral over C is on the segment Rex = ±e with e — >• + , from /j, 2 T/M — /try > Imx > on the first 
Riemann sheet and from fi 2 T/M — fxy < Imx < on the second Riemann sheet. Therefore, we have to add 
to Bq an additional term — 2 i 7r A A with 



AA = A{fi) - A(fi), 



A{t) = [l + i—)t-i 



1 



1/2 



(65) 



Note that in the limit r,7^0we have /3 — > f3 and this additional term vanishes. Furthermore, A(j3) = (3 
and A(f3) — f3 c , with f3 2 = 1 — Am 2 /s p . Therefore, we reproduce the correct result of Eq. (|52"|) . The recipe is, 
therefore, replace In with ln~ in the integrand but deform the integration contour in order to avoid crossing 
of the positive imaginary x-axis when this would occur. 

In summary, our result with a simple example we observe that 



T . Analyt.Cont. T . _ 
Ll n I > Ll n 



L Wi( z ) ¥= 



dx 



Li n (x), 



(66) 



since deformation of the integration contour is required for the general case. 



Im z 



Im z 



Re z 



Re z 



Figure 3: Analytical continuation of a Bo -function as seen in the z -plane with a cut along the positive real axis 
between R- and R+ (Eq. (|60[0 . In the first part the integration path reaches the point 7_ fEq. (|60[) ) and continuation 
after z = I- is in the second Riemann sheet. In the second part, where /_ < R- continuation must be, once again, 
in the second Riemann sheet; therefore the integration path which has moved into the lower half-plane must be 
deformed to cross the cut before moving once more into the lower half-plane (but on the second Riemann sheet). 
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6.5 Narrow width approximation 

The practical implementation for higher point (or higher loop) functions presents a formidable technical 
problem, due to the higher dimension of the a: -space; more will be explained in Section 16.61 but, for this 
reason, we have also considered analytical continuation in narrow- width-approximation (hereafter NWA). 
Here we replace In with In - (or ln + ) at the integrand level and do not perform any deformation of the 
integration (hyper-)contour. The resulting expression is expect to have a range of validity given by T -C M. 
Numerical investigation of the Higgs complex pole shows that NWA returns reliable results when compared 
with the exact expression. The rationale for analytical continuation in NWA is based on the fact that, as 
we are going to show, all higher-point (higher-loop) functions admit integral representations with integrand 
of logarithmic nature (one- loop) or, at most, of poly- logarithmic nature (multi-loop) . 

Consider now the extension to complex variables of an arbitrary scalar three-point function Co (in NWA) , 
defined by 

C = / dxi f 1 dx 2 V~ 1 ~ e/2 (x 1 ,x 2 ), (67) 



Jo Jo 

where n — 4 — e and V is a quadratic form 



V(xi,x 2 ) = axf + bxl + cxix 2 + dxi+ex2 + f- iO = x t Hx + 2K t x + L, (68) 

whose coefficients are related to the internal masses and the external momenta by the relations Hij — 
— Pi ■ Pj,L — m\ and 



Ki = ^(pi-pi + ml-m?), K 2 = ^(P-P-p 1 - Pl +ml-ml), (69) 



with P = pi + p 2 . Let us define the usual Bernstein - Sato - Tkachov (hereafter BST) factors (see Ref. [2"T] ) 
as B3 = L — K l H^ 1 K and BST co-factors X = — H^ 1 K. It is convenient to introduce special notations, 
X = 1, X3 = 0, and V(i i + 1) to denote contractions, i.e. 

V(6l) = V(l, Xl ), V(12)=V(a;i,xi), V(23) = V( Xl ,0). (70) 

In this way we obtain a simple integral representation 



C ° = ~k{l + f dXl if 1 ^2lnF(a; 1 ,x 2 )-i^(X J -X J+1 )lnl/(z r TTl).]}. (71) 



When some or all the invariants are complex, P 2 = — s P with s P = M — iT M and mf = /i| — iji pi (in 
realistic cases, e.g. decay of an unstable particle, p\ 2 are real) we define 



V- = V 



r,7.=o 



(72) 



which includes the — i prescription and write 

Co = \\ + I d Xl [ dx 2 ln~ (V(xi,x 2 ) ; V-(xi,x 2 )) 
m L2 J lj 

1 2 

- 2 £ {Xi ~ Xi+l) ln ~ ( v ( ii + 1 ) ; v-(ii + i))] }. (73) 

i=0 

For instance, with P 2 = —M 2 + i M T, pi 2 = and mi j3 = 0, m 2 — p? — i fi 7 we find that when r/7 < fi/M 
\nV must be continued to the second Riemann sheet for < x 2 < (^7) /(Mr). 

Starting with an integral representation of a three-point function where the integrand is the logarithm 
of a polynomial in parametric space is the safest way of performing analytical continuation; of course, going 
beyond NWA requires contour deformation but even the latter admits a consistent numerical implementation. 
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Nor should one fail to notice 't Hooft and Veltman emphasis, in their seminal work [5^, on this subject: 
they put up warning signs about continuation of their result to complex momenta. 

For higher point functions (L = D,E,F, . . . ) we apply the BST algorithm [30 as many times as it is 
needed to produce logarithms in the integrand and proceed by replacing In with In - , 



L= d{x}hi( X ({x})-iO)^ d{x} In" (74) 

where {x} is the xi, ... ,x n simplex and {x}c(x) ^ s the path that avoids crossing the positive imaginary 
X-axis. NWA amounts to the identification {x}c( x ) = {x}. 

For multi-loop integrals other functions must be extended, e.g. we will use Eq. (|53p , with similar results 
for all generalized Nielsen polylogarithms [26 



'S'n,p( zUST ; z - T ) — 5'n,p( 2lJST ) + 



^ k\ 

k=l 



S n ,p~k{ zUST ) / J Tj Sn—j,p— fc( 2;UST ) 
3=0 3 ' 



0(4 ST -1) 0(^ ST ), (75) 



which is derived by using ln + in the integral representation of the generalized Nielsen polylogarithms. Since 
it has been shown that multi-loop diagrams can be written as integrals of multivariate generalized Nielsen 
polylogarithms [31] our recipe gives the analytical continuation to all orders in perturbation theory, but one 
does have to be careful in one respect: using familiar relations such as splitting of logarithms should be done 
with a grain of salt. 



6.6 Analytical continuation and contour deformation 

Exact analytical continuation at the integrand level can be performed by deforming the integration 
contour into the complex parametric space (for a general treatment see Ref. [28j). In this case we need the 
general definition of In - given in Eq.(|37|). 

To illustrate contour deformation we consider, once again, the case of a Bo -function with equal (complex) 
internal masses. If 

M>2fi nT-M~/>0, (76) 

the function \{x),x € [0, h] crosses the positive imaginary axis (the branch cut of In - ). To avoid crossing 
we deform the x -integration into 

1) * = ijj0t, 

2) x =l t + i^B(l-t), (77) 
with t € [0,1] and 8 a free parameter. For x^ 1 ) we require Imx^H^) < 0,V£ € [0,1]; this is possible if 

B < /3max, With 



r2 



8 = 1 -^- 



1 + A 1 + 4 



M 3 

For we require that Rex^HO = corresponds to Imx' 2 '^) < 0, which requires 8 > /3 m i n , where B, 
is the largest, real, solution of the following equation 



(78) 



For /3 m in < B < /3 max we have that x'" 1 ' 2 '' never cross the positive imaginary axis. Furthermore, we compare 
^(!. 2 )(0) with X*- 1,2 '(r) at fixed t S [0, 1] and replace In — > In - when (always at fixed t) X (F) crosses the 
negative real axis for some value T c . In the example that we are considering, illustrated in Fig. 01 everything 
is particularly simple since Imx(Rex) is always a straight line but our recipe works, as well, in the general 
case and allows for a straightforward algorithmic implementation. 
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Figure 4: Example of contour deformations in computing a scalar two-point functions with equal (complex) internal 
masses and complex p 2 . 



For a general recipe of contour deformation, we proceed by analyzing the case where a Feynman diagram 
can be written as: 



/ dxi ■ ■ ■ dx n V] Ai(xi, 
Jo 



,x n ) lnT^(a;i, . . .,x n ), (80) 



where Vi are multivariate polynomials in x\, . . . ,x n , at most quadratic in each variable (note that all one- 
loop diagrams can be written according to Eq. tjSU]) . see Ref. [37]); actually the procedure works as well when 
each Vi is a quadratic form in, at least, one variable. For each term in the sum, we select one variable x = xi 
(among x\, . . . , x n ) and study the analytical continuation (assuming that Im[V] rG ai masses < 0) 

\nV->\n~V with V = ax 2 + bx + c, (81) 

where a, b, c are polynomials in the remaining Feynman variables. The idea is to deform only the x integration 
contour into the complex plane (when needed) while keeping all other variables (x\, . . . , Xi-\, a^+i, . . . , x n ) 
on the real axis. We define 

a — a r + i a,i b = b r + ibi c = c r + ia x — u + iv. (82) 

The real and imaginary parts of V are then given by: 

ReV^ = a r u 2 — 2 7 aiiiv — a r v 2 + b r u — biV + c r 

= a r {u — u c ) 2 — 2 a% (u — u c ) (v — v c ) — a r (v — v c ) 2 + S r , 

ImV = a,i u 2 + 2, a r u v — a,i v 2 + bi u + b r v + ci 

— a, (u — u c ) 2 + 2 a r (u — u c ) (v — v c ) — a.; (v — v c ) 2 + Si, (83) 

where we introduced the following auxiliary variables: 

a * b r r •* a * b2 

5 = 5 r +i8i = c- —r^. (84) 



2 a 2 ' 4a 
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The curves ReV = and ImV — are hyperbolas with center in x c . We also define an auxiliary function, 



U = —di ReV + a r ImV = 2 \a\ 2 (u — u c ) (v — v c ) + a r 6i — a t 5 r . (85) 

The curve U = is again an hyperbola with center in x c and asymptotes parallel to the u and v axes. The 
branch-cut of In - V in the x complex plane is defined by: 

i. T3 is n P t ^ n \ ^ eV = & C/ > 0, if a r > 0, , Q , 

cut : ReV — & ImV > <J=> < „ TA _ „ TT ^ „ „ (86) 

~ \ ReV = k U < 0, if a r < 0. v ; 

First we study the intersection of the curve ReV — with the real u axis. If A = b 2 r — 4 a r c r < 0, the 
hyperbola ReV = never crosses the real axis and there is no need of contour deformation (this is the first 
case of Fig. [5]). If A > 0, the intersections are given by: 




-br ± \/A 

ReV = Qkv = Q =*> u = uf = — . (87) 

If both solutions are not in [0, 1] a distortion is not needed (second plot in Fig. [S]). Even if Uq G [0, 1], it 
can happen that the intersection occurs for ImV < (as in the third case of Fig. [S]) and there is no need of 
deformation to avoid the cut. In order to understand whether it occurs, we can study where the zeros of V 
{ReV — ImV = 0) are lying. This system of equations has always two and only two solutions x± = u± + iv±, 
whose real and imaginary parts are given by: 

a* 5 



u± =u c ± y/\a\ + oy, v± = v c ± sign (c^) \J\a\ - a r , <r = - -^—rx. (88) 

2 ar 
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Note that these points are also solution of the equation U — and (because of the simple form of U) we can 
conclude that: 

ReV=0&U>0 <S> ReV=Q k\ V ^ V+ li u>Uc ifa r >0, 
cut: ReV=0&ImV>0 e> { \v v_iiu<u c (gg) 

Rel/=0&C/<0 & ReV=0&\ V r V+ f U>Uc ifa r <0. 

I U > U_ II u<u c 

At this point we have all information to fix the new integration contour, starting from x — and ending at 
x = 1 without crossing the branch-cut. Of course, as long as the cut is not crossed, all integration contours 
are equivalent and give the same result: it may fairly be said that we have some freedom in defining the 
deformation and that, at the same time, we can control the correctness of the result by using different paths. 
The general situation is depicted in the fourth plot of Fig. [5] and the new integration contour is defined by 
seven segments: 

-ai t, 

—ai + ift t, 
—ai (1 — t) + a c t + ift, 

a c + * ft (1 -«)+&«, 
a c (l-t) + (l + a 2 )t + i/3 2 , 
l + a 2 +ift(l-*), 

(l + oa)(l-t)+t. (90) 

The coefficients ot\, a 2 , ft and a c can be fixed according to the principle of the minimal deformation to 
avoid crossing the cut. This gives the following conditions: 



(1) 


X 


(2) 


X 


(3) 


X 


(4) 


X 


(5) 


X 


(6) 


X 


(?) 


X 



a\ > 


\u-\ 


if 


< < 1 


& u- < 0, 




Q l 


= 


otherwise, 


ft > 
ft <- 


\v-\ 
\v-\ 


if 


< it (7 < 1 


^ w_ > & a r 
V- < Sz a r 


> o, 
< o, 


ft 


= 


otherwise, 


a 2 > 


Kl 


if 


< u+ < 1 


& u+ > 1, 




Q 2 


= 


otherwise, 


ft < - 

ft > 


K'+l 
Kl 


if 


< u£ < 1 


^ v+ < & a r 
w+ > & a r 


>o, 

<0, 


ft 


= 


otherwise, 


a c = 


u c 




< Uq < u£ < 1, 












a c = 





if 


Uq < < u£ < 1, 






a c 


= 


otherwise. 


a c = 


1 




< Uq < 1 < u+ , 













(91) 



The case where all coefficients vanish corresponds to non-deformation: in this case the paths in Eq. (j90|) are 

(1), (2), (3), (4) x = 0, (5) x = t, (6), (7) x = l, (92) 

and refer to < x < 1 on the real axis. Note that the case A < (where it J are not defined) belongs to 
this class. This general recipe for contour deformation works also in special cases where not all segments 
are needed. For example, in the fifth plot of Fig. [5] the new contour consists of only four segments: using in 
Eq. ([90| the conditions of Eq. (|9"Tj) (a 2 = f3 2 = 0,a c = 1), the last three segments reduce in this case to the 
point x = 1. 

We can now consider in this framework the example of a Bo function with two equal masses. In this case 
we have: 

V(x) = -s P x(l-x)+m 2 = (-M 2 +iMT)x(l - x) + f i 2 -iftj, (93) 
which corresponds to 

1 r a 2 M 2 r MT 

Xc = 2' 4 = fl ~ ~Y ~ 1 r 7 + ~T~ 
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= 1 a*c 1 ^ MM + 7I 1 u /J - M7 
CT 8 2|a| 2 8 2MM 2 + T 2 2MM 2 + T 2 ' 1 J 

Since a r > 0, the cut crosses the segment [0, 1] when <7j < (implying that v_ > 0, < 0), a situation 
which occurs for — M7 > 0. It can be verified by explicit calculation that, in this case, we always have 
< u± < 1, which corresponds to the situation depicted in the last diagram of Fig. [5] and the deformation 
requires five segments {ot\ — a.2 — 0, i.e. the first and the last segment in Eq.([90p reduce to a point). 



6.7 Differential operators in a complex domain 

The procedure of analytical continuation at the basis of Eq. (|80[) deserves an additional comment: how 
can we apply a differential operator at the integrand level, in order to get Eq. (|8H|) ? As an example we 
consider the integral of Eq. (|67[) on which we want to apply the BST algorithm: 



C = 



dx 



dyx 



-l+e/2 



(95) 



where x is a quadratic form. As we know Eq. (j80[) follows from BST functional relation; in order to apply 
the BST algorithm in the complex domain (\ 6 C[x, y}) we introduce 



X 



exp [a In \) 



(96) 



and distort the integration path so that it never crosses the positive imaginary axis of x ( or t ne negative 
imaginary axis of x + )• The (— ) analytical continuation of Co is defined by 



dxdy 



A=0 



x(x,y) 



-l+e/2 



(97) 



where A(a;, y) = is the implicit equation for the integration contour. In practice we change variables, 
x = cti t + fii (t £ [0, 1]) with i = 1, . . . , n, n being the number of segments needed to avoid crossing the 
cut (e.g. see Eq. (f90"]) ). The BST functional relation [27] for quadratic forms and the corresponding linear 
differential operator are 



x"(M)=2>0, N , [d x ]) x (H) , 



V 



where [x] = X\ , . . . ,x n . Consider the following integral, 



F_ = 



dz 



zi r 



1 

B 



X(z) 



1 



2( M 



— ^ (x t - Xi) d Xi 



»=i 



(98) 



(99) 



where T is a curve connecting Zi and zj which never crosses the positive imaginary axis of x f° r z € T. Let 
zj be in the second quadrant and z% outside of it; let zq be the point where Imx(zo) = 0, Rex(^o) < 0. 

Thanks to the — prescription the integrand in Eq.(|99p is a continuous function of z (for z € T) and we 
can write 



F = F W + F ( 2 ) = 



dz 



zi r 



X{z) 



dz 



z r 



X(z) 



(100) 



In the first integral [x]_ — X M and we can apply the BST relation of Eq. (|M|) . In the second one we find 



p(2) 



dz 



z r 
/ 

zo r 



X(z) 



dz exp{- 2 i it {p, + 1)} x M (z) 



z r 



dz exp{~2iTr(fi+l)}V{p,,z,d z ) x^ +1 {z) 



dzV (/i, z, d z ) 



zo r 



X(z) 



(101) 
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where we have used exp{— = exp{— 2in([i + 1)}, showing extension of the BST algorithm into the 
second sheet. Note that the BST relation in the second equality of Eq. (|101[) is of pure algebraic nature and 
that the integration over T will always be parametrized in terms of a T(t) : t £ R so that x G C[t\ and 
VeC[t] <d t >. 

Each segment of Eq. ([90|) is of the type considered in or F_ and therefore we can apply the 

BST algorithm to . It goes without saying that this is the correct procedure, instead of applying BST 
first and performing analytical continuation only in a second step. The result reads as follows: 

n J r&i'x+fti ^ r a ij q 

Cp = ^ a-i Tp Ji^jdxj dy]xTx(aiiX+0i,y)- « / dxA H ln_ Xij + y + A, (102) 

i= 1 JQ J j_ ^ JO " 

Xn = x{oti + j9i,x), Xi2 = xWi,x), Xi3 = x(ot % x + I3 % ,a t x + fa), xa = x{ol%x + fa, 0), (103) 

-Aa = a ' + ' , ^4i2 = - — , ^ 3 =jy — y, A i4: = Y, aa=a>i+Pi, a i2 = f3 i: a i3 = a l4 = I, 

(104) 

where B% is the BST factor and X, Y are the BST co-factors. 

With this example we have shown how to apply differential operators when complex momenta are present: 
first the analytical continuation has to be performed together with the deformation of the integration contour 
and just at the end the differential operator can be correctly applied. 

In conclusion we have shown a practical implementation of the concept that the pole at the mass of 
a stable particle can move into other Riemann sheets where it describes an unstable particle. It is worth 
noting that all cases encountered so far, where both the internal masses and the Mandelstam invariants are 
complex, have never been discussed in the literature, although this step represents the logical extension of 
the complex-mass scheme, allowing for a meaningful introduction of pseudo-observables. 



7 Including QED(QCD) corrections 

In this section we will consider the inclusion of QED(QCD) corrections, both virtual and real. The choice 
s = s H in Eq. (f24| is dictated by the request of a gauge independent definition of pseudo-observables and 
follows, once again, from Nielsen identities. Consider a final state where the inclusion of real QED(QCD) 
corrections is mandatory in order to obtain an infrared finite quantity, e. g. i — > H — > 66. Here, at one- 
loop, we have wave-function renormalization factors for the external fermions and vertex corrections; the 
QED part generates, in the so-called (e , nib) regulator scheme (dimensional regularization for the infrared 
and masses for the collinear limit), a simple infrared pole and double as well as simple collinear logarithm. 
According to our recipe the QED(QCD) vertex correction should be evaluated at complex Higgs momentum 
squared. Let us define 

8a = x H ixl, /3 c 2 = l-4^, /3 2 = l-4^f. (105) 
The residue of the infrared (virtual) pole reads as follows 

P \P 2 , o 1 1 , - # 



i?virt — Re — 



— + 2- — 



ln ~ it+T' (106) 



The infrared pole from real emission originates from the end-point singularity in the phase space integration, 
where P H = pb + and P 2 is arbitrary but real, unless one is willing to extend the phase space definition 
into the complex plane where 6 functions are defined in terms of contour integrals |23) . Using the most 
obvious choice, namely P^ — ~n 2 H we obtain for the (real) infrared residue 

i? rC ai = (/3 2 + l) ln^. (107) 
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Therefore, as expected, cancellation of infrared divergences is spoiled by the need of defining virtual correc- 
tions at a complex value of s. However, the fact that Z H 1 (s) Vf(s) is gauge parameter independent only 
at the complex pole does not exclude a gauge independent sub-set of corrections that can be evaluated at 
arbitrary s. Consider the situation at the one-loop level; here, in front of the O (a) QED corrections we use 



Z H = 1 and the one-loop vertex, F^? ED , is gauge independent for all values of s. If we introduce 



Z H = 1 



9 



16 7T 2 

for the wave-function renormalization factor, we can use 

-,2 



sz H , 



S (H c -> bb) 



9 m b 



2 M„ 



1 



16 ir 2 



V~™{s H )--5Z H (s H ) 



(108) 



(109) 



thus preserving gauge invariance without spoiling infrared safety. Following a well established convention it 
is also convenient to define a deconvoluted pseudo-observable where QED(QCD) corrections are subtracted 
according to theory. 

There is an intriguing alternative for the treatment of QED(QCD) corrections; since the definition of the 
Higgs boson mass is not unique, we could keep it as a free parameter, M H . Then, cancellation of infrared 
poles at the one-loop level requires 



(110) 



and i? v irt = ^rcai- Therefore, there is a value of M H which is infrared safe, 

,2 



Ml = \s H \{ 



1 



■HIT. 



-1/2 



O 



(111) 




Figure 6: Examples of mixed electroweak-QED two-loop diagrams contributing to H — > bb. The solid lines, attached 
to the Higgs boson (dash-line) represent Z/(j>° or W/(f> fields. 

The main question in establishing the consistency of the procedure, definition of QED(QCD) correc- 
tions and their subsequent deconvolution, is about the extendability to higher orders. Using the following 
expansions 

00 2n 00 2n — 1 

n— 1 7i—l w 

and working at O (g 5 ) we will have terms like 

5zPv£ );Q * D , (113) 

which are of the mixed type, electroweak-QED, and where Z H cannot be evaluated at arbitrary values of s. 
In Eq.CEH y_W ;QED is the one-loop QED triangle contributing to H — >• bb. 

However, we also have mixed two-loop diagrams, as given in Fig.[6l There is a well-known identity which 
allows us to extract the infrared behavior of these two-loop diagram, in terms of the product of two one-loop 
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+ IR finite 



Figure 7: Infrared decomposition of a mixed electroweak-QED two-loop diagram. Here the scalar case is presented, 
i.e. the spin structure has been completely neglected and, for instance, the wavy line represents a scalar massless 
line. 

vertices plus an infrared finite reminder (see |32j for the explicit decomposition in the scalar case). The 
decomposition for the scalar case is illustrated in Fig. [7] where the external lines in both one-loop vertices 
are on-shell. By scalar we mean those contributions that do not have powers of the integration momentum 
in the numerator. 

As we have seen, the identity holds at the amplitude level, reflecting the factorization of virtual infrared 
corrections and the fact that virtual infrared poles are always coming from Co -functions (the scalar ones). 
These identities follow from the fact that any diagram with an infrared photon line of momentum q, t + K , 
where K is a certain combination of external momenta as well as of the other loop momenta, gives an 
infrared divergence equivalent to the same diagram evaluated at qi = —K. We thus see that the infrared 
decomposition into products of tensor integrals times infrared Co -functions follows trivially. The explicit 
form of infrared factorization, at the amplitude level, is illustrated in Fig. [8] which shows a class of diagrams 
contributing to the two- loop amplitude for H — > bb. 



v 




+ IR finite 



u 

Figure 8: Infrared decomposition of a mixed electroweak-QED two-loop amplitude. Solid lines represent Z, (j>° or 
W, <j> particles; Dirac spinors for the external lines are included and Cj R is the scalar, infrared divergent, three-point 
function. 

When added to the contribution coming from 8Z^ , the combination 

^™( S )-1<>( S ) (114) 

arises naturally in front of an infrared Co- Therefore, our recipe will be to evaluate Eq. (|114l) at s = s H while 
keeping the remaining QED-like Co functions (the infrared divergent ones) at s = fi 2 H . The difference with 
the original diagram is non resonant and mixes with infrared divergent background contributions, e.g. from 
boxes. 

For one-loop real emission we have diagrams as illustrated in the l.h.s. of Fig. |H] where we use p 2 H = —s H 
and where the infrared singularity arises from the end-point of phase space integration which is controlled 
by a real Higgs boson momentum. The corresponding amplitude is gauge independent by construction, a 
fact that can be easily seen in the infrared divergent soft approximation, the r.h.s. of Fig. [9j where we have 
introduced the eikonal factor 

Jcik(p) = -Qb^-, e(k)-k = 0, (115) 
p ■ k 

e being the photon polarization. The vertex correction, first diagram in the r.h.s. of Fig. [9j when summed 

— 1/2 

with Z H (g) LO gives a gauge invariant contribution if both are evaluated at the Higgs complex pole. 
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Figure 9: Examples of the infrared decomposition of O (gr 3 ) electroweak diagrams with real photon emission. The 
last term in the r.h.s of the equation is the corresponding eikonal factor of Eq. (|115[) . 



There is another example where the introduction of QED corrections seems to be controversial if only 
internal masses are made complesQ- Let us consider the pseudo-observable T {H — > W + W~) 1 with on-shell 
external vector bosons and internal complex masses: the infrared behavior of the one-loop corrections is 
reducible to a scalar vertex 

Co {-Ml , -M 2 w ,-s,s w ,0, s w ) , (116) 

where the difference s w — Ml (s w being the W complex pole) acts as an infrared regulator, removing the 
infrared virtual pole, 

Co (-Ml , -Ml ,- s , Sw ,0,s w )=^- In §^ In ° w - ^ + . . . (117) 

Pw S Pw t S 

where (3 W = 1 — 4 Mil s. If we continue the external masses the result is instead 

Co (s w , -s w ,-s,s w ,0, s w ) = — - — ln~ ^ c "' ^ | \ + IR finite, (118) 

Pew ^ Pew -1 e 

where P2 W = 1 — 4s w /s. Let us consider a realistic example, e.g. gg 4f of Fig. [TJ for the complete 
process there is no problem at all because a photon attached to an internal W boson line cannot be infrared 
divergent. However, the goal is a breakdown of the full process into three components, one of which is the 
pseudo-observable T (H — > W + W~); in order to define T (H — > W + W~ (-f)) it is important to control the 
cancellation between virtual and real infrared corrections and, in this case, the extension to external complex 
masses is more than an option. 



8 Schemes 

For processes which are relevant for the LHC and, in particular, for H — > bb, 77, gg and gg — > H etc, we 
define three different schemes and compare their results. The schemes are: 

• the RMRP scheme which is the usual on-shell scheme where all masses and all Mandelstam invariants 
are real; 

• the CMRP scheme [5] , the complex mass scheme with complex internal W and Z poles (extendable to 
top complex pole) but with real, external, on-shell Higgs, W, Z , etc. legs and with the standard LSZ 
wave-function renormalization; 

^^This point was raised by Thomas Binoth in one of our last conversations. 
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• the CMCP scheme, the (complete) complex mass scheme with complex, external, Higgs (W, Z, etc.) 
where the LSZ procedure is carried out at the Higgs complex pole (on the second Riemann sheet). 

The introduction of three different schemes does not reflect a theoretical uncertainty; only the CMCP scheme 
is fully consistent and comparisons only serve the purpose of quantifying deviations of more familiar schemes 
from the CMCP scheme. 



9 Numerical results 

In this section we examine the numerical impact of computing Higgs pseudo-observables at the Higgs 
complex pole (on the second Riemann sheet of the S -matrix) . We use the parametrization s H — [i H — i fi H j H 
for the Higgs complex pole where now \i H is an input parameter and j H is computed in the standard model. 
In this case \i H plays the role of an input parameter while we prefer ~pr H of Eq. (1181) as conventional definition 
of the Higgs boson mass. The results are compared with on-shell pseudo observables. 

As input parameters for the numerical evaluation we have used the following values 

M w = 80.398 GeV, M z = 91.1876 GeV, m t = 170.9 GeV, T w = 2.093 GeV, 

G F = 1.16637 x 10~ 5 GeV" 2 , a(0) = 1/137.0359911, a s (M z ) = 0.118, T z = 2.4952 GeV. 
For the W, Z complex poles we use 

s v = /i 2 -i/iv7v, /4 = M v - T 2 V , /v 7v = M v T v {^~\ J^j ■ ( n 9) 

In computing H — > gg(gg — > H) we have used a running a s (fi H ) (CMRP) or a s (~p H ) (CMCP). Results for 
the computed "f H are collected in Tab. [TJ 
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0.637 


1.448 


4.373 


39.829 


7h [GeV] 


0.051 
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0.391 
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5.069 


40.847 



Table 1: Standard model prediction for *y H in GeV as a function of fi H ■ The first entry corresponds to a real on-shell 
top quark mass, the second entry to a top quark complex pole derived from T t = r" LO = 1.31 GeV and the last entry 
to a top quark complex pole derived from F t equal to the experimental upper bound of 13.1 GeV. 

For the evaluation of all one-loop functions in the CMCP scheme where, sometimes, a continuation to the 
second Riemann sheet is required we used both the analytical results and the exact numerical integration; 
the two in-house (independent) libraries return results in excellent agreement (typically on the sixth digit 
on one- loop percentage radiative corrections). 

9.1 Numerical differences between the CMRP and CMCP schemes 

For a better understanding of comparisons we define weak corrections to H — > bb as 

A wcak = y/2 (Cp art + B part +R), II 2 - Mh {A + ll) 1/2 , (120) 

separating the corrections into a part coming from three-point functions(C par t), two-point functions(_B part ) 
and a rational term(i?). It is worth noting that there are, in general, strong cancellations among the three 
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contributions: for instance, at 120 GeV we have a C par t of —8.233 % from CMCP exact while the bracket in 
Eq.dU]) is -0.790%. 

Differences between the two schemes are roughly of O (j H / n H ), as expected, and become significant 
above the tt -threshold where the Higgs boson width becomes larger and larger. Since the width of an heavy 
Higgs boson is large it is natural to investigate the goodness of the separation of the production stage from 
the decay process. In general these stages are not independent and may be interconnected by radiative 
effects. Our results confirm the theorem of Ref. [33] : radiative effects are not enhanced in totally inclusive 
pseudo-observables with respect to the naive O {^h/ Hh) argument, unless the Higgs boson is very heavy, in 
which case this ratio is large (at [i H = 500 GeV it reaches 29%) and typical cancellations in the total weak 
correction factor are disturbed, increasing the effect. 

To further understand the differences between the two schemes for high values of /i H , we recall the well- 
known fact that the Higgs wave-function renormalization shows an inverse f3 -behavior at the W, Z threshold. 
In the two schemes, exactly at threshold, we will have 
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(121) 



where B — W, Z . The parameter that regularizes the divergence is therefore 7 B — 7^/2 with some sizable 
effect around the ZZ threshold. To analyze differences between the CMRP and CMCP schemes we fix /i H 
and compute "f H ; then we use Eq. ([120[ ) and compare results with the limit j H = 0. Results are given in 
Tab. [2] where we see variations induced by a finite j H . 





Cpart 


-^part 


R 


tot 



0.03 


-3.673 
-3.760 


-1.999 
-1.990 


+4.514 
+4.009 


-1.658 
-1.741 



0.29 


-0.308 
-0.986 


-0.130 
+0.974 


+3.450 
+2.714 


+3.058 
+2.702 



Table 2: Variations (in percent) at \ih 
according to Eq. (|120[) . 



300, 500 GeV in the components of the total weak corrections to H 



bb 



9.2 Testing the NWA approximation 

In order to analyze the quality of the NWA approximation (Section 16. 5[) we have considered the pure 
one-loop weak corrections to the decay width H — »• bb. It turns out, that up to \i H = 250 GeV (where 
Jh/Hh — 0.011) the approximation is very good, less than 1% (of a ~ 1% correction). Note that analytical 
continuation of three-point functions in the exact CMCP scheme is required above 220 GeV. 

9.3 Complete set of results 

As far as the Higgs boson production cross section in gluon-gluon fusion is concerned we find that the 
effect of replacing the on-shell scheme (for the external Higgs boson) with the complex-pole one is completely 
negligible (around 3 — 4 per mill) for low values of the Higgs mass, a fact that is largely expected. Only 
for higher values, say starting from the tt threshold, where "f H becomes larger and larger, we reach sizable 
differences, above 10% and rapidly increasing. 

• H -> 77, gg, gg ->• H 

A detailed comparison of predictions in the CMRP and CMCP schemes for H — > 77, gg is shown in Fig. [T0l 
The partial decay width H -» gg is shown in Fig. [TT1 where we compare the CMRP(=RMRP) and CMCP 
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Figure 10: Comparison of predictions in the CMRP and CMCP schemes for the decays H 
and for H — > gg (red, dashed line). See Section [8] for the scheme definitions. 



77 (blue, solid line) 




300 320 340 360 380 400 420 440 460 480 500 



IL H [GeV] 

Figure 11: Comparison of the decay width T(H — > gg) in the CMRP (dashed line) and the CMCP (dotted line) 
scheme in the high mass region. The effect of a complex top quark pole in CMCP (with a top total, on-shell, width 
of 13.1 GeV) is given by the blue, solid line. See Section[8]for the scheme definitions. 



schemes. Similarly we compare the partial decay width H — > 77 in the RMRP, CMRP and CMCP schemes 
in Fig. [H 

The relatively large effects in T(H — > 77) or T(H —> gg) at large values of fi H are still compatible 
with the naive 0{^ h /^h) argument. Consider Fig. [TUlfor T(H —> 77); for instance, at fi H = 365 GeV we 
have 7ff/^ ff = 0.065 with a variation between the CMRP and CMCP schemes of 13.3% giving a correction 
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Figure 12: Comparison of the decay width F(H — > 77) in the RMRP (dotted line), the CMRP (dashed line) and 
the CMCP (solid line) scheme. The effect of a complex top quark pole in CMCP (with a top total, on-shell, width 
of 13.1 GeV) is given by the blue, dash-dotted line. The red, dash-double-dotted line corresponds to a top width of 
1.31 GeV. See Section[8]for the scheme definitions. 



factor of 2-f H /fi H . Clearly, the large increase in the Higgs boson width for increasing values of \i H makes it 
questionable to use a perturbative description for the Higgs-resonant part of pp — > X when we have a very 
heavy Higgs boson. 

The relevance of this result is clear: if a light Higgs boson is not discovered, one of the goals of LHC will 
be to exclude a standard model Higgs up to 600 GeV [37] ; already at 500 GeV we have 

0CMCP(S3 H) 4-11 /-,01\ 

— = 1.64, parton level, (122) 

CcMR P (gg -> H) 

comparable to the effect of NLO QCD corrections. We have computed also cr{pp — > H) in the two schemes 
using MSTW 2008 LO parton distribution functions (PDF) [34]. The ratio is given in Fig. fT3| for different 
values of s; in this figure we present 

<Tcmcp(pp ->■ H) ^ 123 ^ 
CTcmrp {pp -> H) ' 

where the numerator is evaluated at ~p H (Eq.(fl"8|)) while the denominator corresponds to [i H . Here we use 

dL" dL99 f 1 dx . . rr \ 
a(pp^ H) = a Q r H — — , — — = / — g(x,fj, F ) g [-,fi F , (124) 
dr H dr J T x \x J 

where t = fi\(jr H )/s, ao is the parton level cross section and g is the gluon PDF (with factorization scale 
fi P = fi H for CMRP and n F = pZ H for CMCP). 

Note that this ratio is only an indicator of the (large) size of the effect since, for a realistic value of the 
cross section, NLO(NNLO) QCD corrections should be included, see Ref. [3S] for updated cross sections at 
the Tevatron and the LHC. 

In Fig. [H] we show the corresponding cross section for ^fs = 3TeV, including an estimate of the uncer- 
tainty induced by varying renormalization and factorization scales (kept equal for simplicity); [i H /2 < fj, R — 
Hp < 2[i H for CMRP and Jl h /2 < fi R = ^ F < 2~p H for CMCP {jl H is given in Ea.(fl~gJ)). In Fig.Elwe do not 
include the uncertainty associated to PDFs. 
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Figure 13: The ratio ocmcp/ccmrp f° r the production cross section pp — > H, as a function of fin, for different 
energies, y/s = 3TeV (red, dotted line), y/s = lOTeV (blue, dashed line) and y/s = 14TeV (black, solid line). The 
cross sections are computed with MSTW2008 LO PDFs with factorization scale fj,p = /j, h for CMRP and fip = Jj-h 
for CMCP. 



The ratio between the two cross sections is stable under these variations, e.g. is between 1.4546 and 
1.4574 at s/s = lOTeV and /j, h — 500 GeV. The production cross sections for y/s = 10, 14 TeV are shown in 
Fig. US] note that here we use a s (M z ) = 0.13934 to be consistent with the LO PDFs [34] . 

In order to better understand the numerical differences in the three schemes we show in Fig. [16] a scalar 
two-point function with two internal Z masses. In the RMRP scheme both the internal masses and the 
Mandelstam invariant /i H are real; in the CMRP scheme we replace Mj with the corresponding complex 
pole, s z ; finally, in the CMCP scheme, also the invariant becomes complex and equal to s H . In Fig. 1161 we 
vary \i H between 100 GeV and 500 GeV where ^ H = 146.89 GeV is huge; here deviations between CMRP and 
CMCP become sizable. This simple example shows the general features of one-loop corrections in the three 
schemes; CMRP - CMCP smoothly interpolate the RMRP results around normal thresholds and, when y? H 
becomes larger and larger, CMCP starts deviating from RMRP - CMRP. 

We also consider the expression for the amplitude A (H — > 77) which can be split into a part containing 
only Co -functions and a rational term. We write 
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12 2 
i- A . 



If we introduce auxiliary variables, 

2 

2 / 2 2 \ 1/2 



2 ' 



—2 ' 
/4 



(125) 



(126) 



the amplitude can be written as A — A c + R where 



A c = --x t — C (0 , , -s H ,m t ,m t , m t ) + 6 (x H - 2 x w ) C (0 , , -s H , s w , s w , s w ) , 

i?= _16^ L+ x £L+6 _ (12?) 

O X w X\Y 

A comparison for the real and imaginary parts in the CMRP and CMCP schemes is shown in Fig. |T7J 
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300 320 340 360 380 400 420 440 460 480 500 

^ H [GeV] 

Figure 14: The production cross section pp — > H at yfs = 3TeV for CMRP (red, wide-dashed line) and CMCP 
(blue, solid line). The shaded areas surrounded by the dashed lines give the scale uncertainty obtained by varying 
/ia/2 < fi a = fip < 2 fi H in the CMRP scheme and JZ H /2 < fi R = fi F < 2~p H (Eq.{T8]l) in the CMCP scheme. We 
have used MSTW2008 LO PDFs. 

• H^bb 

Results for H — > bb arc shown in Tab. [3] for the pure weak percentage one-loop corrections in the three 
schemes and for a wide range of Higgs masses; in Tab. 2] the electroweak (weak + QED) percentage one- loop 
corrections are given. For weak corrections the results for H — > 6& are presented graphically in Fig. [T5] the 
figure shows cusps at the tt -threshold due to the fact that the top quark mass is kept real. The origin of the 
cusps is in a Bq -function with p 2 fixed at the complex Higgs pole. The size of the cusps can be related to the 
large Higgs width at the -threshold as illustrated in Fig. [T9l where we analyze Rc Bq(— s h ; m t ,m t ) around 
the It threshold. Here the solid line corresponds to 7„ =0, whereas dash-lines correspond to increasing 
values of ^ H ■ As one can see the limit ^ H ^ is continuous and there is no artifact due to the analytical 
continuation. The wide-dashed blue line of Fig. [TO] corresponds to a finite value of "f H and to a complex top 
pole (with an on-shell width of 13.1 GeV). 

As it is evident the introduction of a complex top quark pole completely cures the shape of the corrections. 
Therefore, this gives further evidence to using the CMCP scheme, at least from a theoretical point of view 
(the top quark total width is, unfortunately, poorly known). It is worth noting that the numerical effect 
given by the blue curve should be interpreted as an upper bound on the effect of a top quark complex pole 
since the experimental result is an upper bound, T t < 13.1 GeV at 95% C.L. Note that, from theory it follows 
= 1.47 GeV and r~ LO = 1.31 GeV. 

Finally, we observe that the O {^ h /^h) effects, which can reach several percent at large values of fi H , have 
a modest effect on all those processes which start at O (g 2 ) (the effect being on NLO corrections) whereas 
the effect is considerably larger for those processes, e.g. H — > Jj(gg), which start directly at NLO (O (g 6 ))- 

10 Conclusions 

In this paper we have shown how to continue Feynman integrals into the second Riemann sheet, in 
a way that can be easily implemented in any program aimed to compute pseudo-observables related to 
Higgs physics at Tevatron and LHC. Pseudo-observables give, in a natural way, the possibility of translating 
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1— Vs=10TeV 



300 320 340 360 380 400 420 440 460 480 500 

H„ [GeV] 

Figure 15: The production cross section pp — > H at y/s — 10 TeV for CMRP (red, dashed-dotted line) and CMCP 
(blue, dotted line); at y/s = 14 TeV for CMRP (cyan, dashed line) and CMCP (black, solid line). We have used 
MSTW2008 LO PDFs with factorization scale fi F = /i„ for CMRP and ^ F = jZ H for CMCP. 




Re B (arbitrary units) 



Figure 16: The scalar two-point function with two real internal masses (RMRP scheme) Mf or two complex internal 
masses (CMRP - CMCP schemes) s z . B is computed in the RMRP (red dash-line), CMRP (blue dot-line) and 
CMCP (black solid-line) schemes. The Mandelstam invariant is /i| in the RMRP - CMRP schemes and s H in the 
CMCP scheme, with 100 GeV < fi H < 500 GeV. See Section[8]for the scheme definitions. 

experimental data into a language that has a direct connection to unambiguous theoretical calculations. 
Using our framework one can freely compute quantities (otherwise non-existing) like Higgs production cross 
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H H [GeV] 



Figure 17: The Co part of the amplitude for H — > 77 and the corresponding rational term of Eq. (|127[) with 100 GeV < 
/iff < 300 GeV. The black, solid (CMCP) and red, dashed (CMRP) lines give the real part whereas blue, dotted 
(CMRP) and magenta, dash-dotted (CMCP) lines give the imaginary part. The cyan, wide-dashed (orange, dash- 
double-dotted) line gives the rational real(imaginary) part of the amplitude without appreciable differences between 
the schemes. The imaginary part of the rational term is always small and negligible. See Section [8] for the scheme 
definitions. 




150 200 250 300 350 400 



H H [GeV] 

Figure 18: The weak one- loop radiative corrections to H — > bb in the RMRP scheme (red, dotted line), in the CMRP 
scheme (black, dashed line) and in the CMCP scheme (black, solid line). The effect of a complex top quark pole in 
CMCP (with a top total, on-shell, width of 13.1 GeV) is given by the blue wide-dashed line. See Section [8] for the 
scheme definitions. The result corresponding to F t = r" LO = 1.31 GeV has no appreciable difference compared to 
the one at T t = 0. 
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M H [GeV] 


5 W RMRP [%] 


S w CMRP [%] 


S w CMCP [%] 


S w CMCP [%] 




massless 


massless 


massless 


massive 


120 


-0.7890 


-0.7904 


-0.7904 


-0.7837 


130 


-0.9557 


-0.9572 


-0.9573 


-0.9509 


140 


-1.1978 


-1.1986 


-1.1986 


-1.1922 


150 


-1.6215 


-1.6146 


-1.6149 


-1.6078 


160 


-4.2656 


-2.6458 


-2.6690 


-2.6587 


170 


-1.3987 


-1.4914 


-1.4875 


-1.4822 


180 


-2.1989 


-1.9435 


-1.9912 


-1.9858 


190 


-1.0338 


-1.1744 


-1.1590 


-1.1569 


200 


-1.1547 


-1.1967 


-1.1987 


-1.1974 


210 


-1.2452 


-1.2621 


-1.2730 


-1.2723 


220 


-1.3132 


-1.3198 


-1.3379 


-1.3376 


230 


-1.3647 


-1.3665 


-1.3914 


-1.3917 


240 


-1.4047 


-1.4044 


-1.4363 


-1.4370 


OCA 

250 


— 1.4376 


i /i oft: 

— 1.4365 


— 1.4759 


— 1.4769 


260 


-1.4674 


-1.4665 


-1.5138 


-1.5151 


270 


-1.4985 


-1.4981 


-1.5539 


-1.5555 


280 


-1.5357 


-1.5361 


-1.6008 


-1.6026 


290 


-1.5851 


-1.5865 


-1.6604 


-1.6624 


300 


-1.6557 


-1.6582 


-1.7410 


-1.7431 


400 


-0.4736 


-0.4865 


-0.8589 


-0.8633 


450 


+1.4855 


+1.4687 


+ 1.1579 


+ 1.1517 



Table 3: Percentage one-loop pure weak corrections for H — > bb; the first entry is the RMRP scheme, the second 
entry is the CMRP scheme, the third entry is the CMCP scheme while the fourth entry is the CMCP scheme with 
finite m,b (Section |8}. 

section and Higgs partial decay widths. 

An unstable particle cannot belong to the in/out basis of the Hilbert space, nevertheless concepts like 
production or decay of an unstable particle becomes aliases for pseudo-observables that have a well defined 
meaning and a direct relation to measured data. 

In this paper a new scheme is introduced which is the (complete) complex mass scheme with complex, 
external Higgs boson (or, equivalently, any other external unstable particle) where the LSZ procedure is 
carried out at the Higgs complex pole (on the second Riemann sheet). 

Pseudo-observables have been a very useful concept at LEP (e.g. Ref.[36]) and will continue to play an 
important role at LHC, although more difficult if deviations from the SM will emerge; in this case model 
independent approaches are required allowing for the extraction of useful quantities that can be fitted with 
different models. 

The usual objection against moving standard model Higgs pseudo-observables into the second Riemann 
sheet of the S -matrix is that a light Higgs boson, say below 140 GeV, has a very narrow width and the 
effects induced are tiny. Admittedly, it is a well taken point for all practical consequences but one should 
remember that the Higgs boson width rapidly increases after the opening of the WW and ZZ channels and, 
because of this, the on-shell treatment of an external Higgs particle becomes inadequate as a description of 
data if the Higgs is not (very) light. 

Furthermore, most of the experimental plots concerning Higgs physics extend well above, say, 200 GeV 
and, if a light Higgs boson is not discovered, one of the goals of LHC will be to exclude a standard model 
Higgs up to 600 GeV (see Ref. [37] for an exclusion plot of the SM Higgs boson for the various channels as 
well as the combination for masses up to 600 GeV). Already at 500 GeV the ratio CMCP/CMRP for the 
gg H cross section is large and comparable to higher order QCD corrections. 
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M H [GeV] 


<S EW CMRP [%] 


S EW CMCP [%] 


120 


-0.9728 


-0.9729 


130 


-1.1467 


-1.1468 


140 


-1.3941 


-1.3942 


150 


-1.8151 


-1.8154 


160 


-2.8485 


-2.8716 


170 


-1.7039 


-1.7000 


180 


-2.1606 


-2.2083 


190 


-1.3990 


-1.3837 


200 


-1.4262 


-1.4283 


210 


-1.4961 


-1.5071 


220 


-1.5580 


-1.5761 


230 


-1.6087 


-1.6337 


240 


-1.6503 


-1.6824 


250 


— 1.6861 


— 1.7256 


260 


-1.7194 


-1.7669 


270 


-1.7543 


-1.8102 


280 


-1.7953 


-1.8602 


290 


-1.8487 


-1.9228 


300 


-1.9232 


-2.0061 


400 


-0.7754 


-1.1488 


450 


+1.1697 


+0.8571 



Table 4: Percentage one-loop electroweak (weak + QED) corrections for H — > bb with a massive b-quark; the first 
entry is the CMRP scheme, the second entry is the CMCP scheme (Section [8]). 




2 

335 336 337 338 339 340 341 342 343 344 

H H [GeV] 

Figure 19: The Re Bq(—s h ; mt,mt) around the It threshold. The solid line corresponds to a real Higgs mass, 7 H = 0, 
whereas dashed lines correspond to increasing values of ^ H ■ The blue, wide-dashed line corresponds to a finite ^y H 
and to a complex top quark pole (with an on-shell width of 13.1 GeV). 
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On top of all practical implications one should admit that it is hard to sustain a wrong theoretical 
description of experimental data if the correct one is available, independently on the size of the effect. 

Finally, our results show that, above the it -threshold the Higgs- resonant contribution to pp — > X, 
correctly described in the CMCP scheme, is strongly influenced by the large imaginary part of the Higgs 
complex pole and the use of the conventional on-shell description of Higgs pseudo-observables becomes highly 
questionable, even from a numerical point of view. 
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